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ABSTRACT 

We describe a method of estimating the abundance of short-period extrasolar planets based on 
the results of a photometric survey for planetary transits. We apply the method to a 21-night 
survey with the 2.5m Isaac Newton Telescope of ~32000 stars in a ~ 0.5° x 0.5° square 
field including the open cluster NGC 7789. From the colour-magnitude diagram we estimate 
the mass and radius of each star by comparison with the cluster main sequence. We search 
for injected synthetic transits throughout the lightcurve of each star in order to determine 
their recovery rate, and thus calculate the expected number of transit detections and false 
alarms in the survey. We take proper account of the photometric accuracy, time sampling 
of the observations and criteria (signal-to-noise and number of transits) adopted for transit 
detection. Assuming that none of the transit candidates found in the survey will be confirmed 
as real planets, we place conservative upper limits on the abundance of planets as a function 
of planet radius, orbital period and spectral type. 

Key words: methods: statistical - planetary systems - open clusters and associations: indi¬ 
vidual: NGC 7789 


1 INTRODUCTION 


Photometric surveys for transiting extra-solar planets have become 
very popular since the dete ction of the transits exh i bited by the 

S ;-host star HD 209458 iCharbonneau et alfcOOdlBrown et ^ 
. For the first time the radius of an extra-solar planet was 
determined, and the measurement of the orbital inclination lead 
to an estimate of the planetary mass, not just a lower limit. The 
planet HD 209458b was found to have an average density of ~0.38 
g/cm®, significantly less than the average density of Saturn (0.7 
g/cm®), leading to the term “hot Jupiter” for the class of Jupiter 
mass planets with short periods (1-10 d). Transiting planets are 
also very important in that their atmospheric composition may be 


determined from transmission spectroscopy for the brighter host 

stars JCharbonneau et alJ 

2 OO 2 I: iBrown. Libbrecht & Charbonneau 

l2004 IVidal-Madiar et al. 

l2004. Careful modelling 

of the transit 
the presence 
g of the star 

HD 209458b, 
arious groups 

morphology and/or timings may be used to constrain 
of moons or rings and to probe the limb darkenin 
jBrown et alJ200lh. 

Since the discovery of the transiting nature of 
many transit candidates have been put forwards by v 

te.g:IStreet et al J200.3I: iDrake & Cool<l2004l:lBramich et al J200.5^. 

OGLE have been by far the most prolific transit sur 
transit candidates from three observational seasons 

vey with 177 

Udalski et al 

120023: lUdalski et alJ l2002bl; lUdalski et alJ 1200.31: 

Udalski et al 


However, even with the discovery of numerous candi- 
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dates, follow-up observations have confirmed the planetary status 
of o nly six, bringing the total number of transiti ng planets to nine 
(seejBramiche^l 1200.51 and references therein; ISato et al]l2005l: 
iBouchv e'^ni2005 ). This is due to the ubiquity of eclipsing bina¬ 
ries and the many observa tional scenari os involving these systems 
that mimic a transit event iBrow iJ200.3h . Spectroscopic and multi¬ 
band photometric observations are required to rule out the eclips- 
ing binary scenario s and determine the mass of the companion (e.g: 
lAlonso et all2004t) . 

When hunting for new planets, the main advantage of the tran¬ 
sit method over the radial-velocity technique is that many stars may 
be monitored in parallel and to fainter magnitudes, thus probing out 
beyond the Solar neighbourhood. Even though only a small fraction 
(~0.1%) of stars are expected to exhibit a hot Jupiter transit signal, 
by using a large field of view instrument on a crowded star field one 
can monitor enough stars to the precision required to detect a num¬ 
ber of transiting planets. Consequently large charge-coupled device 
(CCD) mosaic cameras are essential to the planet catch potential of 
a transit survey. 

A transit survey produces transit candidates that need follow¬ 
up observations to determine the nature of the transit signals. Can¬ 
didates confirmed as transiting planets add to our database of extra¬ 
solar pla nets and constrain t heir poorly known mass-radius rela¬ 
tionship iBurrows et a ll2004 . To estimate the fraction of stars that 
harbour a planet (the planet fraction) as a function of spectral type 
and planet type we compare the number of transiting planets de¬ 
tected with a calculation of the expected number of transiting planet 
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Table 1. The different subsets of stars used when calculating the expected number of transiting planet detections and false alarms. 


Set Of Stars 

Mass Range 

No. Of Stars 

No. Of Stars With 1999-07 Lightcui've Data 

All Stars 

O.OSMfT^ < Mt < 1.40Mf7i 

32027 

20949 

Late F Stars 

< Mt < l.JOATm 

3129 

2780 

G Stars 

O.SOMq < M* < I.OSMq 

7423 

6711 

K Stars 

O.5OM0 < M* < O.8OM0 

15381 

9690 

M Stars 

O.O8AT0 < M* < O.5OAT0 

6094 

1768 


detections. Even when zero planets are detected (a null result) such 
a calculation places upper limits on the planet fraction. 

In this paper we describe a Monte Carlo method for calculat¬ 
ing detection probabilities (and false alarm rates) of transiting plan¬ 
ets based on photometric data, as functions of various parameters, 
taking into account the following factors: 

(i) Limb darkening effects which tend to make central eclipses 
deeper and grazing eclipses shallower. 

(ii) The effect of orbital inclination on the shape and width of 
the transit lightcurve. 

(iii) The distribution of the photometric data in time and the in¬ 
dividual error bars on each measurement. 

(iv) The signal-to-noise threshold, number of transits, and num¬ 
ber of data points in-transit and out-of-transit required for a detec¬ 
tion. 

We then apply the method to the transit survey described in 
ISramich et aljl200^ (here on referred to as BRA05) to determine 
the expected number of transiting planet detections and place lim¬ 
its on the planet fraction as a function of star and planet type. 

In Section 2 we describe the lightcurve data used in the analy¬ 
sis and in Section 3 we define the detection probabilities and false 
alarm probabilities for an extra-solar planet based on photometric 
data. In Section 4 we present the Monte Carlo method that we used 
to calculate these probabilities and derive limits on the hot Jupiter 
fraction in the field of NGC 7789 as a function of star and planet 
type. In Section 5 we discuss the results and in Section 6 we present 
our conclusions. 


2 THE TRANSIT SURVEY DATA ON NGC 7789 

A transit survey of the field of NGC 7789 was presented in BRA05 
in which ~33000 stars were photometrically monitored in the 
Sloan r' band over three separate runs with dates 1999 June 22-30, 
1999 July 22-31 and 2000 September 10-20. For brevity these runs 
shall be refered to from now on as 1999-06, 1999-07 and 2000-09 
respectively. 

To summarise, in BRA05, Sloan r' — i! colour indices 
were used to construct a colour-magnitude diagram (CMD) and 
thereby identify the cluster main sequence. Fig. Q shows the 
CMD for the stars from chip 4 which was centred on the clus¬ 
ter. Although the cluster main sequence is visible, it is clear that 
most stars in the sample are field stars and not cluster stars. 
A theoretical main sequence model for the stellar mass range 
0.08 Mq < M* < 1.40 Mq was adopted and fitted to the clus¬ 
ter main sequence via magnitude offsets. Using the known cluster 
distance dc = 2337pc and reddening E(B — V) = 0.217, 
and adopting an Einasto la w for the distributio n of the interstellar 
medium in the Milky Way iRobin et alJl200.3|l . a distance d* was 
determined for each star such that the theoretical main sequence 
passes through the star’s position on the CMD. It was argued that 



Figure 1. Instrumental CMD for chip 4 (taken from BRA05) which was 
centred on the cluster. The numbers along the top are masses in units of 
Mq and the numbers along the bottom and right are distances in parsecs. 
Fiducial spectral types are marked on the cluster theoretical main sequence 
for clarity. 


giant stars lie beyond the edge of the galaxy in order to be non- 
saturated in the image data. Hence, it was assumed that each star 
is on the main sequence, and after determining the star’s distance 
d*, the star’s mass M* and radius R, could be read off from its 
position on the theoretical main sequence. 

In this paper we consider the 32027 stars from this data set 
that have a lightcurve from the 2000-09 run, and an assigned dis¬ 
tance, mass and radius. The remaining stars with lightcurves lack a 
colour measurement or were too blue to be assigned a mass and ra¬ 
dius using the adopted theoretical main sequence. We also consider 
the lightcurve data from the 1999-07 run where it exists. BRA05 
searched for transits in the 10-night 1999-07 run and the 11-night 
2000-09 run, 14 months later. The 1999-06 run was too sparsely 
sampled in time to support transit hunting by the adopted search 
technique. 

We are interested in the expected number of transiting planet 
detections (and false alarms) for stars of different masses or, equiv¬ 
alently, spectral types. To facilitate this analysis we consider 4 mu¬ 
tually exclusive subsets of stars which in union make up the set of 
all 32027 stars. These sets are the late F stars, G stars, K stars and 
M stars respectively. TableQshows the number of stars in each set 
and the spectral type/mass ranges to which they correspond. The ta¬ 
ble also includes the number of stars for which 1999-07 lightcurve 
data exists. The m ass ranges for the various spectral types are taken 
from|L an 
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3 DETECTION PROBABILITIES AND FALSE ALARM 
RATES 

In BRA05, a matched filter algorithm was used to search for tran¬ 
sits in the lightcurves by adopting a square “boxcar” shape for the 
transit lightcurve of total width 5At (where At is the transit dura¬ 
tion searched for). This search was based on the transit detection 
statistic: 


q2 

“^tra 


Xconst Xtra 



( 1 ) 


where Xua *^he chi squared of the boxcar transit fit, Xconst 
chi squared of the constant fit, Xout is the chi squared of the boxcar 
transit fit for the Aout out-of-transit data points. Transit candidates 
were chosen using a threshold of 5'tra > 'S'min ~ 10- 

Consider an extra-solar planet of radius Rp, orbital period P 
and orbital inclination i with to as the time of mid-transit. The 
planet orbits a star S, of known mass M* and radius J?*, that has an 
associated lightcurve. We calculate the predicted transit lightcurves 
based on a simple planet-star model: we assume a luminous pri¬ 
mary, linear limb darkening with u — 0.5 and a dark massless 
companion in a circular orbit. Adding this signal into the observed 
lightcurve of the star, we calculate the transit statistic 5tra (Eqn.0 
for each transit event, and then evaluate the following detection 
function: 

A (S, , AjQjji, Rp^ P^ipp) 

j 1 if Stra > for at least predicted transits 
1 0 otherwise 

Using the same procedure as above, but without actually adding the 
predicted transit lightcurve into the observed lightcurve of the star, 
we evaluate the false alarm function: 


^ (S) 'S'ffija, Aj^jjj, Rp, P, i, tp) 

j 1 if Stra > 'S'min for ^s^^t Af^in predicted transits 

1 0 otherwise 

The function D is a trigger function that indicates where the 
data and detection algorithm are capable of detecting a transit of 
the specified type, and F indicates where the data alone suggest 
that such a transit is actually present. The BRA05 lightcurve data 
contains some eclipsing binary stars and possibly transits. Hence 
both D and F are slightly over estimated. 

In the upper panel of Fig. ^ we plot a subsection of the 
lightcurve of star 61377 with an injected 0.02 mag offset of du¬ 
ration 3 h starting at tp = 2451799.5 d. This r' « 18.20 mag 
G star has a mass, radius and distance of 0.96Mq, 0.96J?q and 
3152pc respectively. In the lower panel of Fig. we plot the corre¬ 
sponding periodic functions D (tp) and F {tp) represented by thick 
and thin continuous lines respectively. We adopted S'min = 10’ 
i^min = 1- -Rp = 1-2-Rj, P = 3.338 d and i = 90° for 
this calculation. The function D (tp) attains the value of 1 where 
there is data of sufficient precision to detect a transit. The function 
F (tp) attains the value of 1 where there is data that mimics a tran¬ 
sit signature, and it has clearly been triggered by the injected offset 
in the lightcurve data. 

In Fqn.|^ we integrate out the “nuisance parameters” tp and i 



Figure 2. Upper panel: A subsection of the lightcurve of star 61377 with an 
injected 0.02 mag offset of duration 3 h starting at tp = 2451799.5. Lower 
panel: The periodic functions D {tp) and F (tp) for star 61377 represented 
by thick and thin continuous lines respectively. F (tp) is offset vertically 
from D (tp) by 0.1 for clarity. We adopted Sn,in = 10, An,i„ = 1, 
Up = 1.2Uj,P = 3.338 d and i = 90° for this calculation. 


to get: 


P(det|S,5iT,in,An,jn,Up,P) 


p90 rP 

di dtpf(tp,i)D (S, .S' mi n, An^in, -Rp, F h tp) (4) 
Jp Jp 


where P(det | S, S'pnjpi, Apj^ji^, i?p, P) is the detection probability 
for star S and f{tp,i) is the joint probability distribution function 
(PDF) of tp and i. We assume that the parameters tp and i are in¬ 
dependent, tp is uniformly distributed over 0 < tp < P and 
random orbit orientation for i in the range 0° < i < 90°. Hence 
we may write: 


fdod) 


TT sin i 
180°P 


Combining Fqns.|3and|3we get the detection probability: 


(5) 


P(det|S,5iT,in,An,jn,Pp,P) 

/■ 90 ° fP /'Ksmi\ 

~ J o ^ Yg()°^ J ^ l^mini ^Pi R"i lo) 

( 6 ) 

Using a parallel argument, we obtain an expression for the false 
alarm probability P(fal | S, .Smini Ap^jn, Rp, P) as: 

P{M\S,S^in,Nmin,Rp,P) 

" Jo° Jp'’ ^ 

(7) 


4 MONTE CARLO SIMULATIONS 
4.1 Methodology 

We take the Monte Carlo approac h to evaluating t he detection 
probabilities and false alarm rates iPress et alJll992h . rather than 
attempting to numerically integrate Fqns. and 0 In general, a 
Monte Carlo simulation estimates the probability of an event by 
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Figure 3. (a): Detection probability as a function of 5min for star 61377 with P = 3.338 d and Rp = 1.2Rj. The continuous, dashed and shorter dashed 
lines correspond to = 1, A^min = 2 and Almin = 3 respectively. The false alarm probability as a function of for the same star and planet with 

Almin = 1 is shown by the dotted line. The thick continuous line corresponds to the probability Pt = 0.116 that the planet-star system exhibits transits, (b): 
Detection probability as a function of P for star 61377 with = 10 and Rp = 1.2Rj. Again, the continuous, dashed and shorter dashed lines correspond 
to A^min = 1, Af,^i„ = 2 and A^min = 3 respectively. The false alarm probability is approximately zero for all P and Af,^i„ at this detection threshold. The 
thick continuous curve corresponds to the probability Pt that the planet-star system exhibits transits. 


selecting a large random sample from the parameter space as gov¬ 
erned by the underlying PDF, and then calculating the fraction of 
the sample that satisfy the event criteria. The larger the sample, 
the more accurate the calculated probability. However, the size of 
the sample that may be selected and analysed is usually limited by 
available computing resources. 

For each star S and its corresponding lightcurve, we used 
the Monte Carlo method to calculate P(det | S, -^min’ P) 
and P(fal I S, Rp, H) for a grid in A^rnin’ -^P 

and P. We chose to use a geometric sequence in Rjuin from 
Smin = 2.4toSfnjn = 26.4 with geometric factor 1.025. We also 
chose a geometric sequence in P from P = ldtoP=10d with 
geometric factor 1.004 resulting in 576 period values. We chose 
-^min ^ {l,2,3}andRp = 1.2 Rj. For each grid point we se¬ 
lected a set of = 1000 planets with to and i drawn randomly 
for each planet from the PDF in Eqn.|3 

The grid for R should be fine enough that the difference in pe¬ 
riod A R = Ri+i — Pi between two consecutive grid points 
Pi and Ri+i (where Ri+i > Pi) is such that the difference in 
the number of cycles spanning the duration of the lightcurve is 
less than or equal to a fraction /t of the transit duration (in cy¬ 
cle units). This condition implies that the grid in R should be a 
geometric sequence with geometric factor less than or equal to 
1 + (ftAt/T) where At is the transit duration and T is the du¬ 
ration of the lightcurve. Adopting /t = 0.5, At ~ 2 h for a typical 
transit duration and T = 10.4 d corresponding to the longer 2000- 
09 run yields fiAt/T « 4.0 x 10“®. Hence our choice of grid in R 
is fine enough for the ~35% of stars that have lightcurve data from 
the 2000-09 run alone. For the remaining stars with lightcurve data 
from both runs, adopting the much finer period grid that is required 
makes the Monte Carlo simulations prohibitive in terms of com¬ 
puter processing time. To fully sample the possible period aliases 
introduced by using the adopted grid, the current period value Rcurr 
for each Monte Carlo trial was drawn from a uniform distribution 
on the interval R/\/1.004 < Rcurr < R\/1.004. 

For the 1-10 d planets that we consider, the probability that a 


planet transits is ~10%. Hence, the number of Monte Carlo trials 
that result in a detection is smaller than N-^q by at least a factor 
of 10, leading to a relatively noisy determination of the detection 
probabilities (for Wjyic — 1000, a ~10%). However, when we 
sum the detection probabilities over a large number of stars, typi¬ 
cally ~ 10^ (see Section 4.3), the noise is reduced to an insignif¬ 
icant level (~0.1%). Similarly, it is reduced even further when we 
integrate over various period ranges (see Section 4.4). 


4.2 An Example Simulation 

Let us now consider the complete lightcurve of star 61377. The 
star has 118 data points over 10 nights in its 1999-07 lightcurve 
with a standard deviation of ~0.007 mag and 612 data points over 
11 nights in its 2000-09 lightcurve with a standard deviation of 
~0.010 mag. 

In Fig. 13 we plot the detection probability and false alarm 
probability as functions of the transit statistic detection thresh¬ 
old (Fig. |3(a)t and period (Fig. |3(b)t for star 61377. We used 
NyiQ — 10“ Monte-Carlo trials for this star in order to reduce the 
noise from the simulations for illustrative purposes. The detection 
and false alarm probabilities both decrease with detection thresh¬ 
old Rfnin (Fig. |3(a)) . For this particular star, it can be seen that false 
alarms are very unlikely even at very low detection thresholds. The 
thick continuous line corresponds to the probability Pt = 0.116 
that the planet-star system exhibits transits, calculated from: 


Pt 


Rp + R+ 

a 


0.162 


/ Rp + R» \ 
) 


\Mq) 



( 8 ) 


where a is the orbital radius and the final expression uses Kepler’s 
law. The completeness of the transit search falls rapidly with the 
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Figure 4. (a): Expected number of transiting planet detections for all stars as a function of P for Smin = 10 and Rp = 1.2Rj. The continuous, dashed and 
shorter dashed lines correspond to A^min = 1^ ^min = 2 and A^min = 3 respectively, (b): Expected number of false alarms for all stars as a function of P 
for 5nim = 10 and Pp = 1.2Pj. Again, the continuous, dashed and shorter dashed lines correspond to A^min = 1.-^min = 2 and A^min = 3 respectively. 


adopted detection threshold For = 1, the complete¬ 

ness is ~76.7% at threshold = 4, dropping to ~22.3% for 
^min = 10. 

In Fig. |3(b)| we see the expected P dependence of de¬ 
tection probability on orbital period, but with more detailed period 
structure arising from the detailed time sampling of the observa¬ 
tions. For A^jnjn = 1, orbital periods close to integer values tend 
to have lower detection probabilities since such periods are reso¬ 
nant with the observational gaps during the daytime. Conversely, 
orbital periods close to fractional values tend to have higher detec¬ 
tion probabilities since such periods cover a greater range of orbital 
phases. For example, periods close to ~3.0 d have a detection prob¬ 
ability of ~0.025 whereas periods close to ~2.7 d have a detection 
probability of ~0.036 for this particular star. Fig. |3(b)| also shows 
that as you increase the number of recovered transits required for a 
detection, the detection probability decreases rapidly. 


Rp = 1.2i?j. These quantities have a strong dependence on period 
in the same way as the probabilities from which they are derived 
(see Section 4.2). Fig. |4(b)| clearly shows that, by increasing 
from 1 to 2, the expected number of false alarms is effectively re¬ 
duced to zero for all P. However, introducing this extra constraint 
for a transit detection reduces the expected yield of planets from 
the survey by more than a factor of 2 (Fig. |4(a)^ . Also, setting 
^min ~ 2 is unnecessary since the expected number of false 

alarms for the detection threshold chosen in BRA05 (S'min = 10 
and = 1) is less than 1 for P > 1.34 d, indicating a good 
choice of detection threshold for all except the shortest period hot 
Jupiters. 


4.3 Expected Number Of Transiting Planet Detections 

Assuming that each star S has one planet of radius Rp and pe¬ 
riod P, then the expected number of transiting planet detections 
^det (V, jV min , Rp, P) as a function of star type Y, 

^min’ P is simply the sum of the detection probabilities 

for all stars of the required type: 

^det (V 'S'niini l^mim ^Pi 

= y]P(det|S,5^i„,iV^i„,i?p,P) (9) 

ser 

Similarly, the expected number of false alarms 
A'fal (V, ‘^min’ IV min , Rp, P) is given by: 

iVfal (V, •^min 5 -^min i 5 p) 

= y]P(fal|S,5^i„,iV^in,Ep,P) (10) 

ser 

In Fig.|3 we plot the expected number of transiting planet de¬ 
tections (Fig. 14(3)*} and the expected number of false alarms (Fig. 
|4(b)t for all stars as functions of the period with Sjyijn = 10 and 


4.4 Expected Number Of Transiting Hot Jupiter Detections 

We can now estimate the expected number of detections and false 
alarms for three different planet period ranges, and for different 
star types. We consider the very hot Jupiters with periods of 1-3 d, 
the shorter period hot Jupiters with periods of 3-5 d and the longer 
period hot Jupiters with periods of 5-10 d. Within each period range 
Pi < P < Pm we assume that planets are uniformly distributed 
in log(P). This i s roughly consistent with the results of the radial 
velocity surveys iHeacoxI 19991 : iBramic hl200.5h . 

Assuming that each star S has one planet of radius Rp in 
the specified period range, then the expected number of transiting 
planet detections is obtained by summing over star type Y and in¬ 
tegrating over period P: 

^det (V, 'S'jjjjji, , Rp, Pi, Pm ) 

[Pm / d In P \ 

^ (in (Pm/Pi)) ^ 

( 11 ) 
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All Stars 



(a) 


1 -3 d Very Hot Jupiters 



(b) 


3-5 d Shorter Period Hot Jupiters 



(c) 


5-10 d Longer Period Hot Jupiters 



(d) 


Figure 5. (a): Expected number of transiting hot Jupiter detections (continuous curves) and false alarms (dashed curves) for all stars as functions of S'min with 
= 1 and = 1.2R}. Each cui've is labelled with the hot Jupiter period range to which it corresponds, (b): Expected number of transiting 1-3 d 
hot Jupiter detections as a function of ^min with A^min = 1 ^^d = 1.2R} for various sets of stars. Each curve is labelled with the star type to which it 

corresponds, (c): Expected number of transiting 3-5 d hot Jupiter detections as a function of 5min with A^min = 1 ^nd R^ = 1.2R] for various sets of stars. 
Each curve is labelled with the star type to which it corresponds, (d): Expected number of transiting 5-10 d hot Jupiter detections as a function of ^min with 
i?p = 1.2R} for various sets of stars. Each cuiwe is labelled with the star type to which it coiTesponds. 


Similarly, the expected number of false alarms is given by: 

-^fal (^5 '^min 5 -^min ’ ) 

cPm / d In P \ 

^ Jp, [in {Pm/P i )) ^ 

( 12 ) 

In Fig. |5(a)| we plot the expected number of transiting planet 
detections (continuous curves) and false alarms (dashed curves) for 
all stars as functions of with A'jjjjj, = 1 and Rp — 1.2Rj. 
Each curve is labelled with the period range to which it corre¬ 
sponds. Similarly, each of Figs. |5(b)| - (d) corresponds to a different 
period range in which we plot the expected number of transiting 
planet detections for the F, G, K and M stars in our sample as a 
function of with — 1 and Rp — 1.2Rj. 

For a limited range 7 < < 15, the “curves” for in 

all of Figs. 5(a)-(d) can be approximated by straight lines, allowing 


us to express our results in the form of an approximate empirical 
relationship: 


-Vdet (‘^min) 

A'det (10) 


exp 




(13) 


where A is a constant that may be determined for each set of stars 
V and period range Pi to Pm considered. We note that the transit 
detection statistic Sua is equivalent to the signal-to-noise (S/N) of 
the fitted transit signal and that: 


where Af/fo is the fractional transit depth. Now, since we have 
Stra oc Rp, we may infer that for a fixed S'ti-a the equivalent 
detection threshold varies as oc Rp^. Using this fact, we 
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Table 2. Results of the Monte Carlo simulations for = 10, = 1 and ijp = The number in brackets for and indicates the 

uncertainty on the last decimal place. The upper limits on /p at the significance level a may be scaled with and Rp by using the scaling relation defined 
in Eqn.|21Jtogether with the listed values of A. 


V 

No. Of 
Stars 

Period 

Range 

-^det 

Nfal 

A 

Upper Limit On fp 
Ata = 0.50 

Upper Limit On fp 
Ata = 0.05 

Upper Limit On fp 
Ata = 0.01 

All Stars 

32027 

Id < P < 3d 

325.09(2) 

0.792(1) 

2.18 

0.213% 

0.922% 

1.42% 

All Stars 

32027 

3d < P < 5d 

130.77(2) 

0.267(1) 

2.17 

0.530% 

2.29% 

3.52% 

All Stars 

32027 

5d < P < lOd 

60.59(1) 

0.115(1) 

2.15 

1.14% 

4.94% 

7.60% 

Late F Stars 

3129 

Id < P < 3d 

28.50(1) 

0.223(1) 

4.20 

2.43% 

10.5% 

16.2% 

Late F Stars 

3129 

3d < P < 5d 

10.37(1) 

0.075(1) 

4.10 

6.68% 

28.9% 

44.4% 

Late F Stars 

3129 

5d < P < lOd 

4.38(1) 

0.028(1) 

3.99 

15.8% 

68.4% 

105% 

G Stars 

7423 

Id < P < 3d 

163.57(2) 

0.294(1) 

2.18 

0.424% 

1.83% 

2.82% 

G Stars 

7423 

3d < P < 5d 

65.37(1) 

0.110(1) 

2.20 

1.06% 

4.58% 

7.04% 

G Stars 

7423 

5d < P < lOd 

29.85(1) 

0.051(1) 

2.19 

2.32% 

10.0% 

15.4% 

K Stars 

15381 

Id < P < 3d 

111.76(1) 

0.265(1) 

1.91 

0.620% 

2.68% 

4.12% 

K Stars 

15381 

3d < P < 5d 

46.15(1) 

0.078(1) 

1.94 

1.50% 

6.49% 

9.98% 

K Stars 

15381 

5d < P < lOd 

21.80(1) 

0.035(1) 

1.94 

3.18% 

13.7% 

21.1% 

M Stars 

6094 

Id < P < 3d 

21.26(1) 

0.010(1) 

1.44 

3.26% 

14.1% 

21.7% 

M Stars 

6094 

3d < P < 5d 

8.88(1) 

0.004(1) 

1.51 

7.81% 

33.7% 

51.9% 

M Stars 

6094 

5d < P < lOd 

4.28(1) 

0.002(1) 

1.51 

16.2% 

70.0% 

108% 


may rewrite EQn. ll.3l as: 


(15) 

For Spnin = 10, the range of planetary radii over which this scal¬ 
ing relation is valid is 0.98 i?j < Rp < 1.43i?j. 

The relation in Eqn. 11.51 allows us to scale our results for a 
range of planetary radii that encompasses the radii of the known 
transiting hot Jupiters without having to repeat the Monte Carlo 
simulations. In Table we present the results of the Monte Carlo 
simulations for the various sets of stars and orbital period ranges 
considered, as calculated for = 10, Alpujjj = 1 and 

Rp = 1.2 J?j. We also include the value of A from Eon. llSI for 
each combination of star and planet type. One can see that we ex¬ 
pect to detect ~325 1-3 d very hot Jupiters, ~131 3-5 d hot Jupiters 
and ~61 5-10 d hot Jupiters from the transit survey based on the 
assumption that each star has a single planet of the specified type. 

4.5 Upper Limits On The Hot Jupiter Fraction 

In reality, only a fraction fp of the stars considered in our Monte 
Carlo simulations harbour a planet of a specified type, instead of 
our assumed one planet per star. Hence we must correct our cal¬ 
culations of the expected number of transiting planet detections by 
this factor. We refer to fp as the planet fraction. Since fp is an un¬ 
known quantity that we would like to estimate, we may use the fact 
that the transit survey in BRA05 has most likely produced a null 
result (although this is still to be confirmed) and place a significant 
upper limit on fp. 

First of all we make the assumption that the actual number 
of transiting planet detections X has a Poisson distribution with 
expected value E{X) given by: 

E{X) = /pA^det (16) 

The Poission distribution is defined by: 

p(X = x)= g-E(x) j; £ Ng (17) 

x\ 


(^det (‘5'mini ^p) 

fVdet (10, 


exp 


A 1 - 


10 


E2Rj 

Rn 


For a null result, a: = 0. Using this fact and substituting Eqn. m 
into Ean. llTI we get: 

P(X = 0) = (18) 

In order to obtain an upper limit on fp at the significance level a 
we require: 


P(X = 0) < a 
Hence, from EQns.ll8|and ll9l we derive: 

— In a 


/p ^ 


N, 


det 


(19) 


( 20 ) 


The values of fp that we derive in this manner for a = 0.50, 
a = 0.05 and a = 0.01 are shown in Table|^ For example, for 
a = 0.01, we place an upper limit of 1.42% on the 1-3 d very 
hot Jupiter fraction based on the assumption that such planets have 
a typical radius of 1.2 7?j. Our most conservative upper limits are 
obtained for a = 0.01 and consequently these are the upper limits 
on fp that we consider in the discussion and conclusions. 

Finally, we derive how the upper limit on fp scales with S'fnjn 
and Rp by using EQn. ll5l in Ean. UOl 


/P< 


— In a 

Wdet (10, 1.2i?j) 


A 




5 DISCUSSION 

We have derived relatively stringent upper limits on the abundance 
of hot Jupiters for the field of NGC 7789. The most stringent upper 
limit on fp at a = 0.01 of 0.79% is obtained for 1-3 d very hot 
Jupiters of radius 1.4 Rj. For the Sun-like G stars, we obtain limits 
on fp at a = 0.01 of 2.82% for the 1-3 d very hot Jupiters and 
7.04% for the 3-5 d hot Jupiters (assuming Rp = 1.2Rj). Fig.|^ 
shows plots of the upper limit on the planet fraction fp at a signif¬ 
icance level of 1% as a function of star type (Fig. |6(a)t and orbital 
period (Fig.|6(bjl for an assmued planetary radius of 1.2 Rj. We 





















8 D.M. Bramich and Keith Horne 




(a) 


(b) 


Figure 6. (a): The upper limit on /p at q = 0.01 as a function of star type for the various types of hot Jupiters defined by their period ra nges with 
i?p = 1.2i2j. The dotted line shows the estimate of the 3-5 d hot Jupiter fraction for Solar neighbourhood Sun-like stars from iButler et 'ai] i200(il . (b): The 
upper limit on /p at a = 0.01 as a function of planet type (defined by the period ra nge) with Rr, = 1 .2R\ for various star types. The dotted line shows the 
estimate of the 3-5 d hot Jupiter fraction for Solai' neighbourhood Sun-like stai's from iButler et '^|20o3). 


also show the estimate bv lButler~ 'ai] <2000ll that ~1% of nearby 
Sun-like stars (late F and G dwarfs) host a 3-5 d hot Jupiter, as 
derived from radial velocity observations (dotted line). 

It is interesting to note that although the K stars are the most 
numerous in our star sample (15381 stars), it is the 7423 G stars 
that produce the largest expected number of transiting planets, and 
therefore the strictest upper limits on /p (Tablel^- This is due to the 
fact that the G stars are generally brighter than the K stars in our 
sample, and therefore the corresponding gain in accuracy of the 
photometric measurements outweighs the smaller number of stars 
for which we can search for transits and the smaller transit signal 
for a given planetary radius. 

W e may compare our results directly with those of iButler et alJ 
j200Ch by considering the late F and G stars in our sample and 
the corresponding expected number of transiting planets for the 3- 
5 d hot Jupiters with i?p = 1.2i?j. We expected to detect ~10.4 
such planets around the late F stars in our sample, and ~65.4 such 
planets around the G stars in our sample, an expected total of ~75.7 
3-5 d hot Jupiters. This places an upper limit on /p of ~6.08% at 
the 1 % significance level for these types of star and planet. This is 
consistent with the estimate of the hot Jupiter fraction derived by 
iButler et alJ i200(il) of /p ~ 1% and demonstrates with confidence 
that the hot Jupiter fraction for Sun-like stars in this field may not 
be more than a factor of ~6 times greater than that for the Solar 
neighbourhood. 

By considering the number of planets detected by the radial- 
velocit y technique, and by the transit te c hniqu e from the OGLE 
survey, [Gaudi, Seager & Mallen-Ortielaj i2005h estimate that the 
relative frequency of 1-3 d very hot Jupiters to 3-9 d hot Jupiters 
is O.ISIq J^. Using this result together with /p ~ 1% for hot 
Jupiters, the authors calculate that /p « 0.1% — 0.2% for 1- 
3 d very hot Jupiters. This is consistent with our upper limit on /p 
ata = 0.01 of 1.42% for 1-3 d very hot Jupiters of radius 1.2 i?j. 


6 CONCLUSIONS 

The calculation of the expected number of transiting planet detec¬ 
tions for a transit survey has been discussed in varying levels of 
detail by several authors (e.g:|Gilliland_eLalJ2QQd;|Weldrake_eLaJ 

l2005tlHidas et al 1200.51 : lHoo^^il200^lMochei^^^ll2Q0T ^ 

A requisite for such a calculation is an estimate of the masses and 
radii of the stars in the sample from observational constraints or a 
model for the star population that predicts these properties. 

In this paper, we present a method for calculating in detail the 
detection probabilities (and false alarm rates) of transiting planets 
for photometric time-series data as functions of detection thresh¬ 
old, planetary radius and orbital period. The calculation is based 
on Monte Carlo simulations and requires the properties of the host 
stars to be known, either from observational constraints as in the 
case of stellar clusters or a model for the star population. We have 
shown how to convert these probabilities into an expected number 
of transiting planet detections as a function of star and planet type. 
For a null result in a transit survey, this information can be used to 
determine a significant upper limit on the planet fraction. 

In the case of the transit survey of NGC 7789 presented in 
BRA05 we have derived upper limits on the planet fraction for 
the F, G, K and M stars in the sample and for three relevant pe¬ 
riod ranges of hot Jupiters. We have also derived how these limits 
scale with detection threshold and planetary radius. In BRA05, it 
is estimated that the survey expects to detect ~2 HD 20945 8b-like 
transiting planets or ~4 OGLE-TR-56b-like transiting planets us¬ 
ing simple arguments. Our results indicate that for HD 209458b 
(3-5 d hot Jupiter with J?p ~ 1.4 J?j), and under the assump¬ 
tion that /p « 1%, we also expect to detect ~2 such transiting 
planets. Similarly for OGLE-TR-56b (1-3 d very hot Jupiter with 
Rp ~ 1.2 Rj), and again under the assumption that /p « 1%, 
we expect to detect ~3 such transiting planets. It is encouraging to 
note the agreement between the two methods although the simpler 
method from BRA05 may tend to slightly over estimate the plane¬ 
tary detection rate. We conclude that the transit survey presented in 
BRA05 reached the sensitivity required in order to detect a few hot 
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Jupiters if the abundance of such planets in the field of NGC 7789 
is similar to that of the Solar neighbourhood. 

Improved survey design, mainly by employing a longer survey 
duration, will greatly improve the sensitivity to hot Jupiters. It is 
well known that metal rich star s have a much higher probability 
of hosting an extra-solar planet iSantos. Israelian. & Mavonl2004h 
and hence careful choice of the target star population will increase 
the probability of a detection. Even in the presence of a null result, 
a more sensitive survey will allow the derivation of tighter limits on 
the planet abundance. 
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